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The new perturbation theory for the problem of nonstationary anharmonic os- 
cillator with polynomial nonstationary perturbation is proposed. As a zero order 
approximation the exact wave function of harmonic oscillator with variable fre- 
quency in external field is used. Based on some intrinsic properties of unperturbed 
. wave function the variational-iterational method is proposed, that make it possible 

to correct both the amplitude and the phase of wave function. As an application 
' the first order correction are proposed both for wave function and S'-matrix ele- 

. ments for asymmetric perturbation potential of type V{x,t) = a{T)x^ + P{t)x^. 

I The transition amplitude "ground state - ground state" Wqo{X;p) is analyzed in 

Q^ ' detail depending on perturbation parameter A (including strong coupling region A 

~ 1) and one-dimensional refraction coefficient p. 

Ch ■ 1. Introduction 

Q^, Every successful effort in investigation of some physical problem in strong interaction re- 

gion usually brings new insight in behavior of corresponding system. Thus, the outstand- 
\ ing results, obtained in the solution of one-dimensional model problems in quantum field 

j_j ■ theory, molecular and soUd-state physics [1 — 3] are of great interest. But the investiga- 

^ tion becomes much more complicated when many-dimensional problems are approached. 

In such situation contradictory results are possible even in calculations of the system 
eigenenergy [4] . 

The situation becomes even more dramatic while dealing with nonstationary problems 
with eigenf unctions basis, changing in time. The most important example of it is given 
by model of two-dimensional scattering problem, that is called collinear model for rear- 
rangement [5] . It is one of the simplest and still realistic descriptions of the three-body 
reaction of A + {B, C)n — >■ {A, B)^ + C type. 

This problem can be of principle importance both for the practical applications in the 
theory of chemical reaction for moderate energies and for illuminating the mechanisms 
of rearrangement in multichannel scattering theory. It can be also regarded as a good 
testbed for application of nonperturbative method in realistic few body problem with all 
reaction channels involved. 

In our previous paper [6] it was shown, that in the limit % ^ this problem is ef- 
fectively reduced to that for one-dimensional anharmonic oscillators with the frequency 



Q(t), changing in the field of external force F(t). Numerous quantum-chemical calcula- 
tions of the potential energy surface for most popular three-particle systems have shown 
strong anharmonism in normal coordinate even for principle quantum number as low, as 
n > 3^4. 

Here we propose the method, that seems to be a natural development of the ideas 
of stationary perturbation theory of [3] for the case of nonstationary coefficients of a 
zero-order equation. The method is of variational/iteration type and is essentially backed 
by main ideas of singular perturbation theory - one can expand only such parts of the 
solution, that do not bring singularities, and if one can not find pattern for singular part 
of the solution it is possible to expand it in the exponential. 

Our final goal would be derivation of the scattering matrix representation, that makes 
it possible the simultaneous computation of the probabilities for all pertinent scattering 
mechanisms, taking into account close coupling of many channels, resonance effects and 
chaotic behavior in intermediate state. 

The perturbation theory proposed might be especially interesting for calculation of 
wave packet of scattering system, since harmonic approximation, that exactly propagates 
Gaussian wave packets, is a natural starting point of that theory [7]. 

The suggested development covers also two of the principle questions of the quantum 
theory. First one is connected to the possibility of study of the quantum system behavior 
in the region of close coupling. For that purpose it is necessary to build up a perturbation 
theory, that do not need a knowledge of the spectrum of unperturbed system. Also there 
are several such approaches (see i.e. [3]), they are not generalized to nonstationary case 
and so are not suited for scattering problems. Our method exactly serves for that purpose. 

Second question is closely connected to the limits of quasiclassical description of the 
quantum system evolution. The traditional WKB approach, as it is well known, is limited 
for scattering applications. At the same time there is famous path integral approach that 
makes it possible to describe exactly quantum dynamics in terms of classical trajectories. 
So the basic question is how to improve the WKB approximation to describe all the effects, 
for which traditional approach fails. Our investigations makes it possible to develope the 
representation of the scattering theory via the classical trajectories and so build up the 
intermediate approximate solutions between exact and primitive WKB forms. 

It is well known, that nonstationary form of the equation for quantum oscillator by a 
special transformations can be reduced to stationary form [10], that can be solved exactly. 
In our case unfortunately the use of such approach is not very effective, sine after such 
transformation the perurbation becomes so complicated, that analytic computation of the 
perturbation matrix elements would be unrealistic. 

2. Formulation of the problem 

Let us consider the Schroedinger equation for one- dimensional nonstationary anharmonic 
oscillator 



L * + (x,t)* = 0, 



— OO < X,T < +00, 



(2.1) 



(2.2) 
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with \1/ being the oscillator wave function, A being the dimensionless coupling param- 
eter, and the perturbation potential V{x,t), oscillator frequency ^{t) and the external 
force F{t) having the following asymptotic behavior 

V±ix) = Urn Vix,T), = lirn (](t), lirn F(t) = 0. (2.3) 

Let us assume, that solution of (2.1) satisfies the following boundary and initial con- 
ditions 

lim ^(a;,T) = 0, lim T)(?a;*(x, r) = 0, 

x— >±oo a;— >±oo 

(2.4) 

lim ^(a;, r) = *j„(n; x) exp -i (n + \) ili„r] , 

were ^m(^; ^in,x) is a function of stationary anharmonic oscillator of initial channel. 

The physical meaning of those conditions is obvious - the decay of probability and its 
flow at infinity - and it corresponds to the formulation of the physical problems, discussed 
in Introduction. When perturbation is absent, i.e. with = 0, the equation 

L*o(ri;x,T) =0, (2.5) 
with boundary conditions (2.3), (2.4) has an exact solution (see i.e. [8] ) 

"^foin; X, t) = /o(n; x, r) = K{n] r) exp (ai(T)y + a2{T)y^) Hn{y), (2.6) 
f(n- /7r)i/2\'/' ( } } \ 

y = (^^.n) ICH^ (2.7) 

Lci -livf- + Fv, E{n- t) = ^(n + 1/2), 

with constant being the initial frequency, n being the principle quantum number, 
E{n; r) being the adiabatically changing n-the energy level and the functions C (''")) ^ (''") 
satisfy the equations 

C + fi'(r)C = 0, ^ + (]2(^)^ = F(r), (2.8) 
with following asymptotic and initial conditions 

C+ ('^) = C {t) = ci exp {iO^outr) - C2 exp {-iO^outr) , 

C-(t) = lim C(t) = exp (A„t) , \ci\^ - \c2\^ = 1, (2.9) 

T—> — 00 

Tj (—00) = 77 (—00) = 0. 
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It is well known, that solution of the second equation (2.8) can be constructed on the 
base of solution of corresponding homogeneous equation 

V (r) = [C (r) d* (r) + C* (r) d (r)] , d{r) = ^] dr'C (r) F (r) . (2.10) 

— oo 

Semiclassical type analysis gives the following form of solution of (2.1) 

*(+)(n; X, r) = f{n; x, r) exp (-$(n; x, r)) , (2.11) 

with being the total wavefunction, that is being developed from the n-the asymp- 
totic excited state at r ^ — oo. 

Substitution of (2.11) into (2.1) gives the following equation for unknown functions $ 
and / 

- 1 {d^^f + ^dl^ - XV(x,r)^ f-{Lf- (9.$) (9,/)) = 0. (2.12) 

For further investigation of equation (2.12) it is convenient to expand f{n]x,T) and 
$(n; X, r) functions into power series over A 

oo 

^n;x,T) = Y,X''<^>k{n;x,T), $o(n; r) = 0, (2.13) 

k=0 

oo 

f{n;x,T) = J2X'h{n;x,T). (2.14) 

k=Q 

As in all exponential approximation approaches one term of expansion of $(n;x, r) 
takes into account infinite number of terms of standard perturbation theory. But the more 
we try to take into account in the exponential, the more difficult would be the equations 
for individual members. In our case there is an optimal choice, that makes it possible to 
formulate additional condition for (2.13) , (2.14), and so to determine both functions $ 
and / in (2.11) without ambiguity. Let us start from the equation for "/c"-the correction 

idr^, + ^dl^k + q^\x,T) - fk - id,^k)idJo) + (a;,r)) = 0, (2.15) 

with 

-1 k—l 

= -9 ^ i9x^m){d,^k-m), k>2, 

m=l 

k-1 

qf\x, t) = -Y. [fm{idr^k-m + ;j9^^fe-m - Sk-m,lV- (2.16) 
m=l ^ 

1 k—m—l 

-0 E id,^i)id,^k-m-i)) - {dA){dJk-m)], k>2. 
^ i=\ 

In case oik — \ one has simple relations 
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q['\x,T) = -V{x,t), qr{x,T) = 0. (2.17) 

Note, that the eigenenergies arc not preserved in our problem and so the corrections to 
them in general give no useful information about the system. Thus the problem is reduced 
to the solution of equation (2.15) i.e. to determination of correction to wavefunction of 
anharmonic oscillator. 

3. Construction of the solution 

Before solving the equation (2.15) let us note two important points 

a) The equation is not correctly formulated, as there are two unknown functions in 
one equation, 

b) There is, in general, at least one singular member in this equation. It can be 
explained most easily from analysis of (2.16), taking into account, that /o(n;a;, r) in- 
cludes the Hermitian polynomial. This difficulty is readily solved if the perturbation is of 
polynomial form, as we would suppose in future. 

Such form of perturbation does not seriously influence our discussion, as any smooth 
perturbation can be approximated by a finite part of its Taylor expansion. 

As in every exponential perturbation theory one can solve both above problems by 
choosing additional condition for equation (2.15). Let us start from equations of the first 
order in perturbation. It is convenient to choose 

idr<^i + ^dl<^>i + qi'\n;x,T) + Qi{n;x,T) =0, (3.1) 

Lfi = {d,<^i){dJo) - Qi(n;x,r)/o, (3.2) 

with Qi{n; x, r) being yet unknown function. With additional condition of $i(n; x, r) 
being nonsingular correction, one has that $i(n;a;, r), Qi{n]x,T) and all the terms in 
(3.2) are polynomials of the same order, as qi\n]x,T). By expanding the functions of 
eq. (3.2) in a series of Hermitian polynomials and supposing the absence of the members 
of an order, higher than n, one can correctly determine the coefficients of polynomial 
Qi{n\ X, r). So, it is clear that both above difficulties are thus overcome and the system of 
equations (3.1) , (3.2) is quite correct. The procedure is readily generalized to any order. 

(2) 

As for /c > 2 ql '{n; x, t) ^ 0, so one has 

A;— 1 -j 

Qk{n-x,T) = {{d,^k){dJo)- Lfk- Y.ifm[{idr<^k-m) + 7:dl<^k-m- (3-3) 

m=l ^ 

1 fc— m— 1 

Taking into account, that expression in square brackets is equal to Qfc-ii eq. (3.3) is 
written in the following form 

gfe(n;a:,r) = \{d,^k){d.h)- Lfk- E ifmQk-m - {djmm<^k-m))] U\ (3-4) 

I m=l J 
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and for determination of fk one has the equation 

k—m) fmQk—m) ■ (^■^) 

m=l 

And again if one supposes, that fm, $m and Qm ioi m < k are polynomials, then 
and fk are easily determined and it is possible to find the equation for 

idr^k + ^dl^^k + gi'^ {n; X, r) + Qk{n; x, r) = 0. (3.6) 

Note, that for polynomial Qk{n\ x, r) the solution of (3.6) is also polynomial and only 
high order terms of perturbation potential V{x, r) are important for determination of 
and fm- 

4. Calculation of wavefunction to the first order of perturbation 
theory 

Let us study in more detail the case of asymmetric polynomial perturbation 

V{x, r) = a{T)x' + P{t)x' = ^ ^^(t)^ , (4.1) 

m=0 

with coefficients 6^(r) of the form 

bo{T)=v\f3v + c^), 6i(r) = firn V|C|(4/9r; + 3a), 
h{r) = 3(]-S|Cr(2/5ry + a), hir) = fi-'/'|C|^(4/5ry + a), 

b,ir) = n-:f3\C\\ (4.2) 

and assumption, that q;(t) and (3{t) are slowly varying functions of r. 
Before passing to the solution of (3.1) , (3.2) in accordance to previous discussion let 
us rewrite $i and Qi in the following form 

4 4 

$i(n; x,t) = J2 '"k{n\ T)y^ , Qi{n; x, r) = ^ ak{n; T)y^ . (4.3) 

k=Q k=0 

Thus, the additional condition, introduced for regularization of perturbation approach, 
can be expressed in the form 

4 

/o(n; X, t) ak{T)y^ = Ko((2a2(T)y + ai(T))/o(n; x, r) + 2nK{n; r) x 

k=0 

4 Q 

X K-'{n - 1; r)/o(n - 1; x, r)) ^ ^^^^(t)/-^- L /i, «o = ^:7%- (4.4) 

k=i IsvJI 

Writing down /o(n; x, r) in explicit form and equating on both sides of (4.4) coefficients 
of the polynomials of the same order we may determine the coefficients in polynomial 
Qi{n;x,T) : 
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cro(n; r) = Ko{aiVi + 2nv2 + 2n{n - 1)^4), (Ti{n; r) = 2fi;o(a2f 1 + 01^2 + -^nvs), 

3 2 
(72(n;T) =4^0(02^'+ gait's + ^^^4), c^3(w;t) = 6/«o(a2t'3 + gai'y4), 

(74(77.; r) = 8Koa2'i^4- (4.5) 

Substituting (4.5) in (3.1) one gets a system of nonuniform linear differential equations 
for determination of coefficients in the correction ^i{n]x,T) : 

ivj — Cj{T)vj — dj{T) — 0, Vj — dvj{T)/dT, j — 0, 1, 2, 3, 4. (4-6) 
with functions Cj(r) and dj{T) given by 

C4(t) = AKq, C3(t) = 3^0, C2(r) = 2^0, Ci(t) = Kq, Co(t) = 0, 



d^ir) = 64, c?3(t) = bs, d2{T) = 62 + 2(2n + 3)ko^'4, 



(4.7) 



di{T) = 61 + 3(n + 1)ko?-;3, (io(r) = 60 + (2n + 1)koV2 + 2?7,(n - 1)^0^4. 
Initial conditions for the system (4.5) are 



= 0, = - ido, V] 



J = 1,2,3,4. 



(4.8) 



Note, that " — " parameters correspond to the hmit r — > —00 . It is clear that solution 
of (4.6) must start from j — 4 . In such a way, solution of each equation is presented in 
the following form 



Vj{r) = Gj{T, -00) 



T 

v~ — i j Gj {—00, t') dj (r') dr' 



j = 0,...,4, (4.9) 



with Gj (r, r') being the evolution operator of pertinent homogeneous equation 



G, (r, r') = G-' (r', r) = cxp 



ij Ko{r")dr" 



(4.10) 



Taking into account conditions (4.5) equation (4.4) can be represented in the following 
form 



L fi (n; x,r)-Y^ ij (n; r) fo (n - j; x, r) = 0, 



(4.11) 



with 



64 (n; r) = 2kqV4 



(n-4)! 



1/2 



6*4 (-00, r) , 
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nl 



n 1/2 



_2(n-3)!_ 



G3(-oo,t) , 



62= fi^o [2 (2n - 3) ^4 + V2] 



nl/2 



(n-2)! 



ei (n; r) = [3 {n -l)v3 + 2vi] 



2fn- 1)! 



G'2(-oo,r) , 
1/2 

Gi (-00, r) . 



(4.12) 



In all cases, when n — j -< 0, Cj are equal to zero. The solution of (4.11) is naturally 
represented in such a way 



/i iji] X, r) = ^ Wj (n; r) /o (n - j; r) , Wj {n; r) = Gj (-00, r) u;^ (n; r) . (4.13) 
Substituting (4.13) into (4.11) one gets the linear first-order equations for coefficients 



iwj — Jhq (r) Wj — Cj (n; r) = 0, (n; r) =6^- (n; r) Gj (r, — cxd) , 
with initial conditions 



Cj {n, —00) 



w, — lim w,- (n; r) = 

J r— »— 00 

Solution of this equation is obvious, 



j = l,...,4. 



(4.14) 



(4.15) 



Wj (r) = Gj (r, -00) 



wJ ~ J dr^Gj (—00, r') ej (r') 



(4.16) 



Solutions for higher orders of expansion are constructed in the same way. 



5. Calculation of the transition ^'-matrix for nonstationary an- 
harmonic oscillator 

The fact, that unharmonic oscillator wave functions, determined above, are nonstationary, 
opens new possibilities in view of our previous result [6]. It was shown there, that scat- 
tering in three-body collinear system is effectively reduced to the evolution of anharmonic 
oscillator in external field. And so the above wave functions can be used for computa- 
tion of the scattering matrix. It is not evident how to calculate the scattering matrix 
elements via the oscillator functions. It is possible to show, using the development of [5] , 
that S*— matrix elements for scattering problems can be represented by one-dimensional 
integral in the form, similar to standard representation [9] of nonstationary 5*— matrix: 



Smn= lim {lp}{7n]X,T)llj^ {n]X,T)) , (....) 



00 

/ 



...dx, 



(5.1) 
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with il)f (m; x, r) being the asymptotic wave function of final state. Let us have a look 
at the approximations of exact and asymptotic wave fimctions. For nonstationary wave 
function ip^^^ (n; x, r) with the help of (2.11) to the first order of perturbation theory over 
A one has 



{n; X, t) = [/a (n; x, r) + \fl (n; x, r) 



exp 



1=1 



+ O (A^) , (5.2) 



fx (n; X, t) = exp [-Xvq (r)] /o (n; x, r) , fl (n; x, r) = exp [-Xvq (t)] /i (n; r) . 

Note, that functions /o (n; a:, r) and /i (n; r) are determined by the formulae (2.6) 
(2.7) and (4.13). Expanding the exponential function in (5.2) one has 



ijj^'^^ {n; X, t) = fx {n; x,t) + X 



fx (^; ^> ^) - h (n; X, t) E vi (r) 



1=1 



+ O (A^) . (5.3) 



Now one can use well known formula for Hermitian polynomials xHm {x) = l/2Hm+i 
mHm-i {x) and expand in (5.3) the expressions of type fx {'>^]x,t) into the series over 
the Hermitian polynomials. So one gets 



(n; X, t) = fx (n; x, r) + 



+A 
with 



4 _ 4 _ 

T, wi {n;T) fx{n-l;x,T) - E Up {n;T) fx{n - p; x,t) 

j=l p=-4 



+ O (A^) , 



(5.4) 



u. 



p (n; r) = Gp (r, -oo) (n; r) , (n; r) 



(n — p)! 
2Pn! 



1/2 



Xj,(n;r). (5.5) 



As to Xp (^j 't) functions, they are given by the following expressions 

X-4 (n; r) = ^^;4 (r) , X-3 (n; r) = (r) , 

1 / 3\ 1 3 1 

X-2 (n; ^) = 2 (^^ + 2 j ^ 22^2 ' ^) = 4 + ^) ^3 + ' 

Xo (n; ^) = ^ (^^ + ^ + ^) ^4 (r) + + ^) (r) , (5.6) 
X4 {n; t) ^n{n-l){n- 2) (n - 3) (r) , xs {n; r) ^n{n-l){n- 2) (r) 



X2 {n; r) = 3n (n - 1)^ V4 (r) + n (n - 1) ^2 (r) , Xi (^; 'r) = 2^^'"^ ('^) + ^""i ('^) 
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In the same way wc can get the representations for asymptotic states of anharmonic 
oscillator. So for {out) state one has 



+A 



i/jf (m; X, r) = (p° (m; x, r) + 



E wl (m; r) c/?^ (m - Z; r) - E ("^; r) c/?^ (m - p; r) 



+ 0(A2). 



(5.7) 



Function (m; r) in the above expression is an [out) asymptotic state of harmonic 
oscillator 



tp^j (m; X, t) = ip^j (m; Vt^ut] x) exp 



-i{ni + -] QoutT 



(f^j: (m; fl; x) 



1/2 ^ 

exp ^— -Ox^^ (^VQx^ 



In the above expressions wj (m; r) and (m; r) functions are given by 

w{ (m; r) = w{ [m) exp {-ilQoutr) , 



(n; r) = {n) exp (-ipfio^^r) . 



(5.8) 
(5.9) 

(5.10) 
(5.11) 



As to the coefficients w{ (m) and -u^ [n), they can be obtained at once from (4.16) and 
(5.7) — (5.8) in the limit r —oo after substitution w~ — > Wj and v~ — >■ fj. Now taking 
into account (5.5) and (5.10) and using (5.3) and (5.6) in (5.1) one gets 



Smn (A) — ('S'^n + A + j exp 



with 



4 4 

p=— 4 



1=1 



m(n-l) \"') ^m{n-p)i 

p=— 4 



(5.12) 

(5.13) 

(5.14) 
(5.15) 



Now, starting from (5.12) — (5.15) , it is not difficult to get the analytic expressions 
for transition probabilities for anharmonic oscillator, 



Wmn (A) = exp 



-2Aw 



(+) 



1 + 2Ai?e Q 



WL + 0{X'), (5.16) 



12 
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mn I mn 



lim Wmn (A) . 

A — >0 



(5.17) 
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In the expressions above, the value is a matrix clement of S'-matrix for transitions 
in harmonic oscillator with variable frequency Q (r) in the external field. The matrix 
element in (5.15) is calculated via the generating function [6] (see also [10]). So we can 
give here only the final result 

WL = [^Y' \Hrnniyi,y2)\'exp[-u{l - ^cos29)]. (5.18) 
Here Hmn {vi, Z/2) is the Hermitian polynomial of two variables [10] with 



y,^^u(l-p)e'', y, ^ (^e-'' - ^e'') , 9 ^ ^(5^ + 82) - (3. (5.19) 

z/, p, 5j and (3 parameters are determined from solution of the classical problem for 
harmonic oscillator (2.3) , (2.9) and are given by following expressions: 




^out (1 — p)^^^ 




o^ = e''\l^-^2, C2 = e-^,li^[-^\ , (5.20) 



P = 



d — lim d (r) — \/ve 



T— »+00 



Note, that in the expression proposed both for 5"- matrix (5.1) and for transition prob- 
ability Wjy^n (5.16) we used only the first members of the expansion of exponentials in 
Taylor series over coordinate y. Such approximations for wave function and transition 
probability are effectively used for many interesting applications. In some cases, when 
anharmonic perturbation substantially change spectrum of the problem and thus our 
expansion bases, the leading terms of perturbation must be taken into account. 

It is clear, that for scattering problems unharmonic perturbation is very important for 
higher excited states, especially for rearrangement procceses. At the same time in some 
situations for rearrangement processes even for ground state the effect of perturbation 
can be substantial. To show it let us discuss one particular case. 



6. Application to "ground state-ground state" transition 

We shall demonstrate an application of the method proposed to most simple situation - the 
perturbation of parametric harmonic oscillator by symmetric potential /3(t) ^ with (r) 
being adiabatically changing functions with boundary values /3 {—06) = and /3 {00) = 
P'^ 7^ 0. Note, that transition probabihty for unperturbed oscillator can be obtained 
from (5.16) with 1/ — 0, or by Taylor expansion of generating function of corresponding 
,S-matrix [10] 

were generating function I is equal 
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I {z,, Z2; p) = (1 - pf' exp |1 [vp (zl - zl) + 2yr^ziZ2] } , (6.3) 

with n< = min (m,n) , n> = max (m,n), and {x) being the associated Legendre 
polynomial. 

For many reasons one of most important parameters of the problem is Woo (A; p), that 
measures the probability of change of initial ground state to final ground state. So we 
shall discuss it in some details. From (6.1) — (6.3) one gets 



W^oo (A; p) = exp 



—\v, 



(+)■ 



[1 - 2ReK {p)\ Wo?) (p) , 



(6.4) 



k=0 



with matrix elements S^^^ and S^'j^ given by (6.2) — (6.3), 



(0) 



(6.5) 



c(0) 



1/4 



o(0) 

'-'20 



o(0) 
~'-'02 



V2! 



■\/pSqq , 



c(0) _ q(0) _ 1 ^q{0) 
'-'40 — '-'04 — ^^^^ 00 ■ 



p^r. (6.6) 



For determination of ut2k ^'^d mi'2)j coefficients, that are important for A (p) depen- 
dence and so the dependence of scattering amplitude, one must find coefficients Vj and 

Vj^^ — Vj (00). After integration of (4.9) by parts and taking into account adiabatic 
dependence P (r) one gets 



vj 0, — + -G, (r; -00) / dr'Gj (-00, r') j = 1, 2, 3, 4. (6.7) 



j/i^o (r) 



Now the coefficients Vj^^ are determined from (6.7) taking into account (4.7) — (4.8) 
via integration by parts and averaging over fast oscillations 

4*'^^, (6.8) 

ko= J^m/(|ci|^+ |C2|^). 

In the same way, for coefficients Vj one has 

C4 ^ ^ 'out C2 ^ ^ 'out 

j^out. 

Now using (5.5) — (5.7) from (6.8) — (6.9) one gets the expressions for coefficients u'^2k 
and u'L2k form 

/ 3 /'/3+ ^ + 3 /?+ /l + pV ,R,m 
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The formula (6.4) with the results above, gives the final expression for Wqq 



Woo 

with new notations 



(A;p) ~ - A [l - v+ (p)] (l - ^p)} exp 



(6.11) 



,+ 



+ 



,(+) 



1 + p 



A 



8QI 



-A. 



out 



As it was shown in [10] , parameter p, that measures the excitation of classical oscil- 
lator, corresponds to quantum mechanical refection coefficient of the particle with mo- 
mentum k{x) = ^ [x] . That makes it possible to use well known results from quantum 
mechanics for p. 

As can be seen from fig.l anharmonic oscillator unlike the harmonic parametric one 
in the limit p — > have a transition probability non equal to one. More than that as it is 
seen from (6.11) the dependence of transition probability over A is regular, that makes it 
possible to use the proposed formula up to the values of A of the order of unity. 



Fig. 1. The dependence of "ground state-ground state" transition probability over 
reflection coefficient p and compliny constant A. 

7. Conclusions 

Many important problems in theoretical and mathematical physics are reduced to the 
solution of the equation of nonstationary unharmonic oscillator with different sets of 
[in) and [puC] states. The use of standard approaches to such problem meets two basic 
difficulties: 
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a) Opposite to the stationary situation we do not have the fixed basic set for perturbed 

wave function; 

b) In nonstationary situation the dimcnsionlcss perturbation parameter changes with 
time and can become not small in strong coupling region. 

In this communication, as a generalization of nonlinearization method [3] , we propose 
the way to overcome the above difficulties. The perturbation theory is constructed on 
exact wave functions for quantum harmonic oscillator with variable frequency in external 
field as basis state. This makes it possible, due to some unique intrinsic properties of 
those solutions, to work out the system of two linear equations, that determine the first 
corrections both to the amplitude and to the phase of total wave function. It should be 
noted also, that n-th order correction is also determined by two independent differential 
equations, that are obtained after n iterations. 

Very important computational property of our approach is connected to the fact, 
that any order correction is constructed on the bases of finite number of wave vectors. 
That makes it possible the analytical calculation of the corrections even to the transition 
operator. 

The possibility of simultaneous correction of the phase and amplitude of the wave 
functions, as in a case of stationary problem, gives a regular method of investigation of 
strong coupling region, where perturbation is strong and small parameter is absent. 

The approach is well suited for use in numerous application, of which we point out 
only some most important: 

* corrections to path integral representation of the propagator in field and scattering 
theories, 

* account of kinetics in some problems of solid state and condensed matter physics, 

* anharmonic corrections for scattering matrix in coUinear model of rearrangement 
collisions, 

* anharmonic corrections for the propagator in the wave packet dynamics approach 
to molecular scattering. 

We shall discuss those problems in detail in our future publications. 

Note, that analysis of proposed expressions show, that account for antisymmetric 
terms in potential expansion can cause the nonadiabatic behavior and so nonanalytic 
dependence of Wqq (A; p) over A. Depending on the form of perturbation "ground state - 
ground state" transition probability might have in such a case several maxima. 

Authors gratefully acknowledge the helpful discussions with Dr. Yu. Gorbachev. 
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